Spectral function and conductivity in the normal state of the cuprates: a spin 

fluctuation study. 
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We study the spectral function Ak{u}) and the optical conductivity ai{ij) for a system of fermions 
interacting with their collective spin fluctuations (the spin-fermion model), and apply the results 
to optimally doped cuprates in the normal state. We show that there is no qualitative distinction 
between hot and cold regions in the Brillouin zone - in both cases, there exists a wide range of 
frequencies where the width of the peak in Ak{u!) scales linearly with uj. We demonstrate that cti (uj) 
is inversely linear in uj up to very high frequencies. We argue that these results agree quantitatively 
with the photoemission and optical data for Bi2212. 



PACS numbers:71.10.Ca,74.20.Fg,74.25.-q 

The observed discrepancies between the normal state 
properties of the cuprates and the predictions of Fermi 
liquid theory continue to attract a lot of attention from 
the condensed-matter community. According to a Fermi 
liquid theory, the behavior of all fermionic systems at 
sufficiently low energies is universal and is governed by 
the fact that the quasiparticle damping near the Fermi 
surface scales as w^ or T^ (whichever is larger). This 
should give rise to the w^ (T^) behavior of the width of 
the photoemission peak, to 1/w^ behavior of the optical 
conductivity (Ti(a;) at T = 0, and to the T^ behavior of 
the resistivity p{T). 

Traces of Fermi liquid behavior have been observed in 
strongly overdoped cuprates. For smaller dopings, how- 
ever, the deviations from the Fermi liquid behavior be- 
come substantial, and even at optimal doping, the sys- 
tem behavior in all experimentally accessable frequency 
ranges is qualitatively different from that in a Fermi liq- 
uid pi. Amazingly enough, this non- Fermi liquid behav- 
ior can, to good accuracy, be described by simple linear 
functions of frequency and temperature. In particular, 
the width of the ARPES peak scales with tu for a wide 
range of frequencies, the conductivity is inversely propor- 
tional to w, and the resistivity is linear in T. To which 
extent this behavior survives in underdoped cuprates is 
unclear as below optimal doping, the normal state be- 
havior is masked by the development of the pseudogap. 

There are currently three qualitatively different phe- 
nomenological scenarios for the observed linearities at 
optimal doping. The first is a marginal Fermi liquid sce- 
nario Ig]. It assumes that near optimal doping, there ex- 
ists a quantum critical point of unknown origin, at which 
T.'j^(lu) is linear in uj and independent of k. The linearity 
of the resistivity and inverse optical conductivity then 
follow from a conventional Drude theory. 

The second is a cold spot scenario B. It assumes that 
the quasiparticle lifetime is very anisotropic along the 
Fermi surface and preserves a Fermi liquid form only in 
a narrow range near the zone diagonal where the quasi- 
particle lifetime is the longest. The width of this range 



is assumed to scale with to (or T). Elementary manipu- 
lations show that the Fermi liquid region then yields the 
linear behavior of resistivity and inverse conductivity. 

The third is a magnetic hot spot scenario which re- 
lates the linear behavior of the resistivity with the non- 
Fermi liquid behavior due to strong spin-fluctuation scat- 
tering. This scenario is questionable on general grounds 
as the computation of conductivity implies averaging of 
the quasiparticle lifetime over the Fermi surface |j] , but 
still may be applicable by numerical reasons. 

In the present communication, we adopt a non- 
phenomenological approach and compute the ARPES 
lineshape, optical conductivity and resistivity for the 
spin-fermion model which describes fermions interacting 
with their own collective spin degrees of freedom. Two 
of us argued in a series of recent publications J5J that the 
strong coupling limit of this model captures the physics 
of the cuprates. We will use a previously obtained ex- 
pression for the fermionic self energy S as an input, and 
compute the photoemission intensity for various fc, opti- 
cal conductivity and resistivity. 

The results we obtain partly agree and partly dis- 
agree with each of the three phenomenological scenari- 
ous above. In agreement with the marginal Fermi liquid 
picture, we found that S"(w) is linear in lu and in T for 
a wide range of frequencies everywhere on the Fermi sur- 
face. This linearity in S" causes the linear behavior of 
both the optical conductivity and the resistivity. In con- 
tradiction with the marginal Fermi liquid scenario the lin- 
ear behavior of S"(w) in the spin-fermion model is not as- 
sociated with a closeness to an unknown phase transition 
at optimal doping. Rather, it emerges as an intermedi- 
ate asymptotics in the crossover regime between the two 
physically motivated limits: a Fermi liquid at the lowest 
frequencies, where E" ex w^ and a high frequency limit 
where the system is in the magnetic quantum-critical 
regime, and S" ex ^/uj H. The crossover between the 
two limits is governed by a single parameter - a momen- 
tum dependent typical spin relaxation frequency LOsf{k). 
We find, however, that the crossover region is strikingly 



wide and covers the whole frequency range probed in the 
experiments. 

In agreement with the cold spot scenario, we found 
that the amplitude of the scattering rate is anisotropic 
over the Fermi surface, and is smallest near the zone di- 
agonals. Acoordingly, (i) the slope of S(w) is the smallest 
and (ii) the crossover frequency ujsf{k) is the largest for 
k = kdiag- However, we found that ujsfikdiag) is very low, 
(< 50meV), such that for w and/or ttT probed in ARPES 
and conductivity experiments, the system is outside the 
Fermi liquid regime even for k = kdiag- 

Finally, in agreement with the hot spot theory, we find 
that the damping rate is indeed the strongest near the 
hot spots, and that there is a substantial, although not 
dominant, contribution to the conductivity from the hot 
regions. 

We now turn to the calculations. The fermionic self- 
energy in the spin-fermion model has been obtained be- 
fore ||5|] and is the starting point for our studies in this 
paper. We have 
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where cOm and tUn are fermionic Matsubara frequen- 
cies. Here A ^ ^ is the dimensionless effective cou- 
pling, k = |k — k/is| is the momentum deviation 
from a hot spot along the Fermi surface, and ujsf — 
{3sm(po/ {16tt))vp£^~^ / X, where vp is the Fermi veloc- 
ity, 00 is the angle between the Fermi velocities at k and 
k + Q, and ^ - is the magnetic correlation length. 

The physical meaning of tOsf can be understood by an- 
alyzing Eq. (pi) at T = and k ~ k^s- In this limit, the 
frequency summation and transformation to a real axis 
can be performed exactly and yields Sfc^g (^) = 2 Xu;/(1 + 
\J\ — i\Lo\lLOsj). A simple manipulation then shows that 
LOsf is a crossover scale between a Fermi liquid behav- 
ior at Lj ^ ujsf where I]fc^^(a;) « X(lo + iLi;|cj|/(4Li;s/)), 
and the quantum-critical (^ — oo) behavior at cj ^ ujaf 
where I]fe^^(w) w isigna;(|w| w)^'^, and Q — 4A^Ws/ is 
independent on ^. 

In Fig. |l|a,b we present plots of Efc^^(w,T). We com- 
puted YI'{lo^T) by using the spectral representation for 
(0), and then used a Kramers-Kronig transform to ob- 
tain E'(a;,T). We see that at a given T , S"(w) is linear 
in Lo at intermediate frequencies. At T = 0, the lin- 
ear regime extends between 0.5ws/ and Sw^/. This lin- 
ear behavior has been observed in numerical studies [p|. 
In Fig.nlc we see that S'^' (T) at fixed, small w is lin- 
ear in T above V^QK. Theoretically, the linearity in T 
can be understood as coming primarily from the scatter- 
ing by thermal, classical spin fluctuations (i.e., from the 
n = m term in (0)). In Fig.nld-f we present the results for 
the photoemission intensity /fe(w) = Ak{Ltj)np(Lu), where 
A,iu) = {lM^'l{co)/{ic. -e, + T.'^{u)Y + {Y."{u)f) is 
the quasiparticle spectral function. Fig. Id shows /fe(^) 
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FIG. 1. Theoretical results for the fermionic self-energy, 
Eq. hi and the photoemission intensity. For definiteness, we 
used A = 1, ujsf = 20meV and vp = 3eVA (this yields 
^ ~ 9A). Figs, a and b - T.'{uj) and T,"{uj) respectively 
at various T, indicated on the figures. Fig. c - E"(r) at 
uj = 0. Fig. d - the photoemission intensity vs et- Fig. e- 
the FWHM of the photoemission peak vs oj at e^ = (ECD 
curve). Inset - the FWHM of the peak in fc— space at oj = vs 
T (MDC curve). Fig. f - the velocity of the EDC dispersion 
Vp — vf/{1 + dE'/da;) vs oj. 



for various ej,. Obviously, the width of the peak increases 
with increasing e^. In Fig.np, we plot the full width at 
half maximum (FWHAI) of the peak vs w at Cfe = 
(EDC curve). We see that the width scales linearly with 
frequency over a wide frequency range. In the inset to 
this figure we plot the FWHM of the peak in fc— space at 
UJ — and varying temperatures (MDC curve). Again, 
the width clearly scales linearly with T over a substantial 
T— range. In Fig. If, we plot the velocity measured in the 
EDC dispersion Vp ~ wf/(1 + dYi'^{w) / du) vs frequency. 
Ai UJ = Q, v*p — Vp / {1 + X) . At larger frequencies, the 
renormalization factor decreases due to the fiattening of 
S', and v*p gradually approaches vp. 

We now discuss how these results are modified away 
from the hot spots. For this we notice that at finite 
k — k — khs , Eq. (P can be reduced to the same form as 
at a hot spot, if we introduce the fc dependent coupling 
Xk and spin fluctuation frequency Usf (fc) via 

Afc = A/(l + (fcC)')l/^ c^,/(fc)=c^,/(l-f (fc^^) (2) 

We see therefore that away from the hot spots, the 
effective coupling gets smaller, and the crossover fre- 
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FIG. 2. A comparison between the theory and the ARPES 

data. a,b. A fit of Argonne data |7| for k ~ khs (a) and 

k = kdiag (b); c. Fits of the FWHM of the EDO curves 

to the theory. The circles with error bars are data from [H] 

taken at 115 K and the triangles from lid] taken at 90 K. The 

two theoretical curves are for different coupling constants (see 

text) 



quency ujsf{k) increases. Still, however, at frequen- 
cies/temperatures which exceed ujsf{k), Sj!(w,T) is lin- 
ear in both frequency and temperature. Obviously then, 
for intermediate us, changing the momentum along the 
Fermi surface only affects the overall slopes of S'^'(w) vs 
Lu and T. Furthermore, at optimal doing C ~ (1 ~ 2)a, 
where a ~ 3^ is a Cu — Cu distance, and |fco| ~ 0.37r/a. 
Accordingly, A^ changes between A and (0.7 — 0.4)A be- 
tween khs and kdiag, i-e., the change in the coupling con- 
stant is not dramatic. 

This analysis is indeed an approximate one, and care 
must be taken when applying it to momenta about 
ki? — \i.diag where 0o = '^, and Eq. (n|) is, strictly 
speaking, unapplicable. This complication certainly af- 
fects the fermionic self-energy at the smallest frequencies. 
However, we verified numerically that for frequencies 
uj > LUsf (kdiag), the dominant contribution to EJ,' . (uj) 
comes from fermions away from diagonal, for which (J)q 
can still be approximated by a constant. We therefoe will 
use Eq. (|l|) for momenta both near the hot spots and 
along the zone diagonals with the understanding that in 
the latter case, the result has to be modified at the lowest 
frequencies. 

In Fig. 0, we compare our theoretical results with the 
ARPES data for optimally doped Bi22l2. In Figs.|a and 
b, we fit ARPES data taken by the Argonne group 



to our formulas for k ^ kdiag and k = khs- The val- 
ues Xk^i^^ = 1, uJsfikdiag) = 20meV, and Afe^^ = 2, 
^sf{khs) — 5meV. respectively were used. The value 
of Afe^^ agrees with our previous study of the position of 
the resonance peak at optimal doping Q . To account for 
the background, we added a momentum and frequency 
independent damping 7 = 70meV to ^,"(10). This damp- 
ing likely comes from impurity scattering, but this is not 
well understood and requires further study. 

In Fig. 0c, we compare the theoretical slope of the 
FWHM of the EDC curve with the Argonne data for the 
deviations from the Fermi surface along the zone diag- 
onal. Again we added a constant damping 7 = 70meV 
to E"((x'). We see that the agreement is quite impressive. 
The second line in Fig. ||c is our fit to the Brookhaven |1^] 
and Stanford |l^ data, which also yields a linear in fre- 
quency FWHM, but with a different slope. To fit this 
data, we just have to use A^^.^ « 2.5. This value of 
the coupling is also quite reasonable. A slightly diffent 7 
(55 meV) was used for the second line. We also verified 
that the results for the FWHM are largely insensitive to 
changes in Usf which for Fig. gc is 20 meV. 

ARPES lineshapes at khs and kdiag have previously 
been fitted using the marginal Fermi liquid phenomenol- 
ogy by Sj!(w) = Alo + Bk with /c— independent prefactor 
A |12|. Our analysis shows that such data can equally 
well be fitted by E'^'(a;) = AkOJ + B with a moderately 
A:— dependent A^ and constant B. We emphasize, how- 
ever, that the theoretical reasoning is qualitatively dif- 
ferent in the two approaches. 

We now turn to the conrputation of the conductivity 
(j{uj) = iH(w)/a', where H((u;) is a current-current corre- 
lator. Diagramatically, H((a;) is given by a particle-hole 
bubble with dek/dk in the vertices. In our calculations, 
the typical momenta are comparable to kp, and hence 
vertices can be treated as constants. Also, since the self- 
energy, Eq. (|l|), depends predominantly on frequency, 
vertex corrections to the particle-hole bubble (related to 
dYj/dk by the Ward identity) do not change the physics 
and can be safely neglected. For k near the hot spots, the 
neglect of vertex corrections can also be justified by the 
argument that the velocities at k and k + Q are nearly 
orthogonal, and hence the transport lifetime is the same 
as the conventional lifetime. 

Substituting the expressions for the self-energy into the 
particle-hole bubble and expanding, as before, to linear 
order in deviations from the Fermi surface, one can ex- 
plicitly perform the integration over ek and obtain 
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+ Sfe(w„ -t- UJ„i) ~ i;fc(w„) 



(3) 



where the momentum integration is along the Fermi sur- 
face, I]fe(w„) is given by (|l|), 9(x) is a Theta-function, 
and LUpi is the plasma frequency. We used ujpi ^^ 1.2 x 
lO^CTO^^, similar to that in il3l. 
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FIG. 3. Optical conductivity at optimal doping. For def- 
initeness we used A = 2, Usf = 5meV a. The theoretical re- 
sults for a\(ijj). The experimental data is from [hsl is shown 
in the inset; b. af^ vs lj for T = 907^ and T = ZmK. The 
lines are theoretical curves, and the points are the experi- 
mental data; c. a^^{ui) in the extended frequency range. The 
experimental data is from Ref |13]; e)the resistivity piT) vs 
T. 



The fc— dependence of Sfe(a;) emerges through the mo- 
mentum dependence of Afc and u>sf{k). As this depen- 
dence is not dramatic at optimal doping, it should not 
substantially affect the frequency dependence of n('^). 
To make the computations more transparent, we replace 
Afc and ujsf{k) by averaged, fc-independent values. 

We computed 11" (w) by using a spectral representa- 
tion, and used a Kramers-Kronig transformation to ob- 
tain n'(cj). In Fig. ya we present our theoretical results 
for aiiuj) = Rea{uj) at T=90 and 300 K. The inset of 
this figure contains experimental data from Ref ||l^. In 
Fig. gb, we plot l/ai{uj) for frequencies above lOOOcm"-^ 
and compare it with the data. We see that the theoreti- 
cal l/cri(a;) is linear in frequency. This fully agrees with 
the data. Furthermore, we found that A = 2, compatible 
to that which we used for the spectral function, yields 
a perfect fit to the measured slopes of l/ai{iv) both at 
90K and 300K. We consider this agreement a strong 
indication that spin-fluctuation mechanism captures the 
essential physics of the cuprates. To match the values of 
the conductivity, we, however, again have to add a fre- 
quency independent (but temperature dependent) con- 
stant to l/(Ti. As roughly, l/ai{u!) oc T,"{uj), this proce- 
dure is qualitatively similar (although not quite equiva- 
lent) to adding a constant to T,"{uj). Similar procedure 
has been used in Ref. Q. In Fig. ||c we plot our l/ai{uj) 
in the extended frequency range, up to 8000cto~^. We 
see that the theoretical linear behavior of l/ci extends 
up to very large frequencies of 7000cm^^, where S"(w) 
already curves down from a linear behavior. This exten- 
sion of the linear regime to very high frequencies is also 
consistent with the experimental data of p^ . Finally, in 



Fig. ||e, we present the result for the resistivity p{T). We 
see that p{T) is linear in T above lOQK. This is indeed 
a consequence of the linearity of E"(T, w = 0). 

To summarize, in this paper we used the spin- 
fluctuation approach to calculate the photoemission in- 
tensity, optical conductivity and resistivity in the normal 
state of the cuprates near optimal doping. We found that 
the fermionic self-energy is linear in both lo and T in a 
wide range of frequencies and temperatures. This gives 
rise to a linear frequency dependence of the inverse con- 
ductivity and to a linear temperature dependence of the 
resistivity. We performed quantitative comparisons with 
the experimental data and found near perfect matches of 
the slopes for both the ARPES linewidth and the conduc- 
tivity. We view the results as a strong indication that the 
interaction between low-energy fermions and their spin 
collective degrees of freedom is the dominant scattering 
mechanism in the cuprates. The open issue is whether 
the spin-fluctuation scenario is capable of explaining the 
differences between the temperature dependence of the 
diagonal and Hall conductivities |g,|lj,|l^. This study is 
currently under way. 

It is our pleasure to thank D. Basov, G. Blumberg, 
J.C. Campuzano P. Johnson, M. Norman, J. Schmalian, 
O. Tschernishev and V. Yakovenko for useful discussions. 
We are also thankful to D. Basov for sending us the 
unpublished experimental data. The research was sup- 
ported by NSF DMR-9979749 (Ar. A and A. Ch.) and 
by NSF DMR-96-32527 (R.H.). 



[1] 
[2] 



[3] 
[4] 
[5] 

[6] 

[7] 

[8] 

it 



P.W. Anderson, The Theory of Superconductivity in the 
High Tc Cuprates (Princeton University Press, Princeton, 
New Jersey, 1997). 

C. Varma et at, Phys. Rev. Lett. 63, 1996 (1989), ibid 
64, 497 (1990). CM. Varma, Phys. Rev. B 55, 14554 
(1997). A different but somewhat similar in spirit is the 
scenario based on the assumed closeness to a charge den- 
sity wave instability: C. Castellani, C. DiCastro and M. 
Grilli, Phys. Rev. Lett. 75, 4650 (1995). 
L. B. loffe and A. J. Millis, Phys. Rev. B, 58, 11631 
(1998). 

R. Hlubina and T.M. Rice, Phys. Rev. B, 52, 13043 
(1995). 

A. Chubukov, Europhys. Lett. 44, 655 (1997); Ar. 
Abanov and A. Chubukov, Phys. Rev. Lett., 83, 1652 
(1999). Ar. Abanov, A. Chubukov, and A. Finkel'stein, 



cond-mat/9911445 



T.P. Devereaux and A. P. Kampf, Int. J. Mod. Phys. B 

11, 2093 (1997). 

A. Kaminski et aL, Phys. Rev. Lett. 84 , 1788 (2000). 

A. Kaminski et al, cond-mat/0004482 , 

Ar. Abanov, A. Chubukov, and J. Schmalian, cond 



mat/0005163 



[10] 

[11] 

[12] 
[13] 



[14] 



T. Valla et ai, Science, 285, 2210 (1999). 
Z-X Shen et ai, Science, 280, 259 (1998). 
E. Abrahams and C. Varma, Proc. Natl. Acad, of Sci- 
ences 97, 5714 (2000). 

A. Puchkov, D. Basov and T. Timusk, J. Phys.: Cond. 
Matter 8, 10049 (1996). D. Basov, private communica- 
tion. 

B.P. Stojkovic and D. Pines 
(1997); ibid 5 6, 11931 (1997) 



Phys. Rev. B 55, 8576 
see also A. Rosch, cond- 



mat/9910432 



[15] A. T. Zheleznyak, V. M. Yakovenko, and H. D. Drew, 
Phys. Rev. B 59, 207 (1999). 



